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ABSTRACT 


Through a combination of deep wide-held imaging 
and multi-object spectroscopy, the Next Generation 
Space Telescope will be able to chart with unprece¬ 
dented accuracy the evolution of cosmic structures af¬ 
ter the ‘dark ages’ (z < 5), when galaxies are thought 
to assemble and form the bulk of their stars. In this 
talk I will discuss how NGST observations might de¬ 
termine the history of the stellar birthrate in the uni¬ 
verse and reveal the character of cosmological ioniz¬ 
ing sources at high redshifts. 


1. INTRODUCTION 


The remarkable progress in our understanding of 
faint galaxy data made possible by the combina¬ 
tion of HST deep imaging (Williams et al. 1996) and 
ground-based spectroscopy (Lilly et al. 1996; Ellis et 
al. 1996; Steidel et al. 1996b), has recently permit¬ 
ted to shed some light on the evolution of the stel¬ 
lar birthrate in the universe, to identify the epoch 
1 < z < 2 where most of the optical extragalactic 
background light was produced, and to set impor¬ 
tant contraints on galaxy formation scenarios (e.g., 
Madau et al. 1998c; Steidel et al. 1998). While one of 
the topical questions concerning our understanding of 
the emission history of the universe is the nature and 
redsliift distribution of the recently discovered popu¬ 
lation of sub-mm sources (Hughes et al. 1998; Barger 
et al. 1998; Lilly, this volume), one could also imag¬ 
ine the existence of a large population of faint galax¬ 
ies still undetected at high-,?, as the color-selected 
ground-based and Hubble Deep Field (HDF) samples 
include only the brightest star-forming young ob¬ 
jects. In hierarchical clustering scenarios, a popu¬ 
lation of high-z dwarfs (i.e. an early generation of 
stars in dark matter halos with circular velocities 
Vcirc ~ 50 kms -1 ) is actually expected to be one of 
the main source of UV photons and heavy elements 
at early epochs (e.g. Miralda-Escude & Rees 1998; 
Loeb, this volume). 

While it is hard to predict in detail what the scien¬ 
tific focus of NGST will be ten years from now in a 
rapidly evolving area like physical cosmology, I be¬ 
lieve it is a useful exercise to describe a couple of 
illustrative science programs which well highlight the 
strategic capabilities of NGST in providing a com¬ 
prehensive inventory of the buildup of metals, stars, 


and light in the universe: (a) the determination of 
the rate of Type la and Type II supernova (SN) ex¬ 
plosions as a function of cosmic time from intermedi¬ 
ate to high redshifts; and (b) the study of candidate 
sources of photoionization at early times. Through¬ 
out this paper I will adopt an Einstein-de Sitter uni¬ 
verse (qo = 0.5) with H 0 = 50/i5o kms -1 Mpc -1 . 


2. SUPERNOVAE 


The evolution of the SN rate with redshift contains 
unique information on the star formation history of 
the universe, the initial mass function (IMF) of stars, 
and the nature of the binary companion in Type la 
events. All are essential ingredients for understand¬ 
ing galaxy formation, cosmic chemical evolution, and 
the mechanisms which determined the efficiency of 
the conversion of gas into stars in galaxies at various 
epochs. While the frequency of “core-collapse su¬ 
pernovae”, SN II and possibly SN Ib/c, which have 
short-lived progenitors, is essentially related, for a 
given IMF, to the instantaneous stellar birthrate of 
massive stars, Type la SNe - which are believed to re¬ 
sult from the thermonuclear disruption of C-0 white 
dwarfs in binary systems - follow a slower evolution¬ 
ary clock, and can then be used as a probe of the 
past history of star formation in galaxies (e.g., Yun- 
gelson & Livio 1998). The recent detection of Type 
la SNe at cosmological distances (Garnavich et al. 
1998; Perlmutter et al. 1998) allows for the first time 
a detailed comparison at z < 1 between the SN rates 
self-consistently predicted by stellar evolution models 
that reproduce the optical spectrophotometric prop¬ 
erties of field galaxies, and the observed values. Ac¬ 
curate measurements at all redshifts of the frequen¬ 
cies of Type II(+Ib/c) and la SNe could be used 
as an independent test for the star formation and 
heavy element enrichment history of the universe, 
and significantly improve our understanding of the 
intrinsic nature and age of the populations involved 
in the SN explosions (Sadat et al. 1998; Madau et 
al. 1998a). A determination of the amount of star 
formation at early epochs may be of crucial impor¬ 
tance, as the two competing scenarios for galaxy for¬ 
mation, monolithic collapse - where spheroidal sys¬ 
tems formed early and rapidly, experiencing a bright 
starburst phase at high-z (Eggen et al. 1962; Tins¬ 
ley & Gunn 1976) - and hierarchical clustering - 
where ellipticals form continuosly by the merger of 
disk/bulge systems (White & Frenk 1991; Kauffmann 
et al. 1993) and most galaxies never experience star 
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Figure 1. Evolution of the observed comoving luminosity density at rest-frame wavelengths of 0.15 (dotted line), 0.28 
(solid line), 0.44 (short-dashed line), 1.0 (long-dashed line), and 2.2 (dot-dashed line) /im. The data points with error 
bars are taken from Lilly et al. (1996) (filled dots), Connolly et al. (1997) (empty squares), Madau et al. (1996) (filled 
squares), Ellis et al. (1996) (empty triangles), Treyer et al. (1997) (filled triangle), and Gardner et al. (1997) (empty 
dot), (a) This model assumes a Salpeter IMF, SMC-type dust in a foreground screen, and a universal A 1500 = 0.8 mag. 
(b) This model - designed to mimick a monolithic collapse scenario - assumes a Salpeter IMF and a dust opacity which 
increases rapidly with redshift, ^1500 = 0.09(1 + z) 2 ' 2 mag. 


formation rates in excess of a few solar masses per 
year (Baugh et al. 1998) - appear to make rather dif¬ 
ferent predictions in this regard. By detecting Type 
II SNe at liigh-z, NGST should provide an important 
test for distinguishing between different scenarios of 
galaxy formation. 


2.1. Cosmic Star Formation History 

The emission history of field galaxies at ultraviolet, 
optical, and near-infrared wavelengths can be mod¬ 
eled by tracing the evolution with cosmic time of their 
luminosity density, 

L u y(L„, z)dL u = T(2 + a)?/*L*, (1) 

where y(L v ,z) is the best-fit Schechter luminosity 
function in each redshift bin. The integrated light ra¬ 
diated per unit volume from the entire galaxy popula¬ 
tion is an average over cosmic time of the stochastic, 
possibly short-lived star formation episodes of indi¬ 
vidual galaxies, and follows a relatively simple depen¬ 
dence on redshift. A stellar evolution model, defined 
by a time-dependent star formation rate per unit vol¬ 
ume, a universal IMF, cf>(in ), and some amount 

of dust reddening, can actually reproduce the optical 
data reasonably well (Madau et al. 1998c). In such a 
system, the luminosity density at time t is given by 
the convolution integral 

pu{t) = Peso f l v {t')il>(t - t')dt ', (2) 

J 0 

where l v (t') is the specific luminosity radiated per 
unit initial mass by a generation of stars with age 
t ', and p esc is a time-independent term equal to the 
fraction of emitted photons which are not absorbed 
by dust. The function if(t) is derived from the ob¬ 
served UV luminosity density, and is then used as 



input to the population synthesis code of Bruzual 
& Chariot (1998). Figure la shows the model pre¬ 
dictions for the evolution of p v for a Salpeter func¬ 
tion, A 1500 = 0.8 mag with SMC-type dust, and a 
star formation history which traces the rise, peak, 
and sharp drop of the UV emissivity. For simplic¬ 
ity, the metallicity was fixed to solar values and the 
IMF truncated at 0.1 and 125 M 0 . The data points 
show the observed luminosity density in six broad 
passbands centered around 0.15, 0.20, 0.28, 0.44, 1.0, 
and 2.2 /im. The model is able to account for the en¬ 
tire background light recorded in the galaxy counts 
down to the very faint magnitude levels probed by 
the HDF, and produces visible mass-to-light ratios 
at the present epoch which are consistent with the 
values observed in nearby galaxies of various morpho¬ 
logical types. The bulk (> 60% by mass) of the stars 
present today formed relatively recently (z < 1-5), 
consistently with the expectations from a broad class 
of hierarchical clustering cosmologies (Baugh et al. 
1998), and in good agreement with the low level of 
metal enrichment observed at high redshifts in the 
damped Lya systems (Pettini et al. 1997). 

One of the biggest uncertainties in our understanding 
of the evolution of luminous matter in the universe 
is represented by the poorly constrained amount of 
starlight that was absorbed by dust and reradiated in 
the far-IR at early epochs. Figure 16 shows the model 
predictions for a monolithic collapse scenario, where 
half of the present-day stars were formed at z > 2.5 
and were enshrouded by an increasing amount of 
dust, A 1500 = 0.09(1 + z) 2 ' 2 mag. This results in 
a correction to the rate of star formation by a factor 
of ~ 5 at z = 3 and « 15 at z = 4. The two models 
in Figure 1 produce a significant fraction (~ 50%) of 
the IR background detected by COBE (Dwek et al. 
1998; Fixsen et al. 1998). 
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Figure 2. Predicted Type la and II(+Ib/c) rest-frame frequencies as a function of redshift. The rates are normalized to 
the emitted blue luminosity density. Solid line: SN II rate. Dashed-dotted line: SN la rate with r = 0.3 Gyr. Dotted 
line: SN la rate with T = 1 Gyr. Dashed line: SN la rate with t = 3 Gyr. The data points with error bars have been 
derived from the measurements of Cappellaro et al. (1997), Tammann et al. (1994), Evans et al. (1989), and Pain et al. 
(1997), and have been weighted according to the local blue luminosity function by spectral type of Heyl et al. (1997). (a) 
Model predictions for the merging scenario of Figure la. (b) Same for the monolithic collapse scenario of Figure lb. 


2.2. Type la and II(+Ib/c) Supernova Rates 


Single stars with mass > 8 M 0 evolve rapidly (< 50 
Myr) through all phases of central nuclear burning, 
ending their life as Type II SNe with different char¬ 
acteristics depending on the progenitor mass. For a 
Salpeter IMF the core-collapse supernova rate can be 
related to the stellar birthrate according to 

f„ 125 dmMm) 

SNR n (t) = m r ? 25j (3) 

J 0 1 dmmcpym) 


The specific evolutionary history leading to a Type 
la event remains instead an unsettled question. SN 
la are believed to result from the explosion of C-0 
white dwarfs (WDs) triggered by the accretion of ma¬ 
terial from a companion, the nature of which is still 
unknown (see Ruiz-Lapuente et al. 1997 for a recent 
review). In a double degenerate (DD) system, for ex¬ 
ample, such elusive companion is another WD: the 
exploding WD reaches the Chandrasekhar limit and 
carbon ignition occurs at its center. In the single 
degenerate (SD) model instead, the companion is a 
nondegenerate, evolved star that fills its Roche lobe 
and pours hydrogen or helium onto the WD (Iben 
& Tutukov 1984). While in the latter the clock for 
the explosion is set by the lifetime of the primary 
star, and, e.g., by how long it takes to the compan¬ 
ion to evolve and fill its Roche lobe, in the former it 
is controlled by the lifetime of the primary star and 
by the time it takes to shorten the separation of the 
two WDs as a result of gravitational wave radiation. 
The evolution of the rate depends then, among other 
things, on the unknown mass distribution of the sec¬ 
ondary binary components in the SD model, or on 
the distribution of the initial separations of the two 
WDs in the DD model. 

To shed light into the identification issue and, in par¬ 
ticular, on the clock-mechanism for the explosion of 


Type la’s, it is useful to parametrize the rate of Type 
la’s in terms of a characteristic explosion timescale, 
r - which defines an explosion probability per WD 
assumed to be independent of time - and an explo¬ 
sion efficiency, rj. The former accounts for the time 
it takes in the various models to go from a newly born 
(primary) WD to the SN explosion itself: a spread 
of “delay” times results from the combination of a 
variety of initial conditions, such as the mass ratio of 
the binary system, the distribution of initial separa¬ 
tions, the influence of metallicity on the mass trans¬ 
fer rate and accretion efficiency, etc. The latter sim¬ 
ply accounts for the fraction of stars in binary sys¬ 
tems that, because of unfavorable initial conditions, 
will never undergo a SN la explosion. Possible pro¬ 
genitors are all systems in which the primary star 
has an initial mass higher than m m ; n = 3 Mg (final 
mass > 0.72 Mg, Weidemann 1987) and lower than 
m. max = 8 Mg: stars less massive than 3 M 0 will not 
produce a catastrophic event even if the companion 
has comparable mass, while stars more massive than 
8 Mg will undergo core collapse, generating a Type II 
explosion. With these assumptions the rate of Type 
la events at any one time will be given by the sum 
of the explosions of all the binary WDs produced in 
the past that have not had the time to explode yet, 
i.e. 


SNRi a (t) = 


V fo V>( t')dt ' exp (-T)<j>(m)dm 

t f (f>(m)dm 


( 4 ) 


where m c = max[m m i n , m(t')\, m(t') = (10 Gyr/f') 0 - 4 
is the minimum mass of a star that reaches the WD 
phase at time t', t m = 10 Gyr/ to 2 5 is the standard 
lifetime of a star of mass to (all stellar masses are 
expressed in solar units), and T = (t — t' — t m )/T. 
For a fixed initial mass to, the frequency of Type la 
events peaks at an epoch that reflects an “effective” 
delay At ps r + t m from stellar birth. A prompter 
(smaller r) explosion results in a higher SN la rate 
at early epochs. 
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Figure 3. Predicted cumulative number of Type la and II(+Ib/c) SNe above a given redshift z in a 
Solid line: Type II’s. Dashed-dotted line: Type la’s with r = 0.3 Gyr. Dotted line: Type la’s with 
line: Type la’s with r = 3 Gyr. The effect of dust extinction on the detectability of SNe has not 
models, (a) Model predictions for the merging scenario of Figure la. (b) Same for the monolithic 
Figure lb. 


4' x 4' NGST field, 
r = 1 Gyr. Dashed 
been included in the 
collapse scenario of 


In Figure 2 the predicted Type la and II(+Ib/c) rest- 
frame frequencies are shown as a function of redshift. 
Expressed in SNu (one SNu corresponding to 1 SN 
per 100 years per 10 1 o Lb q ), the Type II rate is ba¬ 
sically proportional to the ratio between the UV and 
blue galaxy luminosity densities, and is therefore in¬ 
dependent of cosmology. Unlike the SN frequency 
per unit volume, which will trace the evolution of 
the stellar birthrate, the frequency of Type II events 
per unit blue luminosity is a monotonic increasing 
function of redshift, and depends only weakly on the 
assumed star formation history. The Type la rates 
plotted in the figure assume characteristic “delay” 
timescales after the collapse of the primary star to a 
WD equal to r = 0.3,1 and 3 Gyr, which virtually 
encompass all relevant possibilities. The SN la explo¬ 
sion efficiency was left as an adjustable parameter to 
reproduce the observed ratio of SN II to SN la explo¬ 
sion rates in the local Universe (SNRn/SNRi a ~ 3.5), 
5% < 77 < 10% for the adopted models. It appears 
that observational determinations of the SN la rate 
at 2 ~ 1 can unambiguously identify the appropriate 
delay time. In particular, we estimate that measur¬ 
ing the frequency of Type la events at both 2 ~ 0.5 
and z ~ 1 with an error of 20% or lower would al¬ 
low one to determine this timescale to within about 
30%. This kind of observations are by no means pro¬ 
hibitive, and these goals could be achieved within a 
couple of years. In fact, ongoing searches for high -2 
SNe are currently able to discover and study about a 
dozen new events per observing session in the redshift 
range 0.4-1.0, and the observations are carried out at 
a rate of about four sessions a year. Since determin¬ 
ing the frequency of SN la with a 20% uncertainty 
requires statistics on more than 25 objects per red- 
shft bin, it is clear that, barring systematic biases, 
those rates will soon be know with high accuracy. 
Also note how, relative to the merging scenario, the 
monolithic collapse model predicts Type la rates (in 
SNU) that are, in the r = 0.3 Gyr case, a factor of 
1.6 and 4.9 higher at 2 = 2 and 4, respectively, with 
even larger factors found in the case of longer delays. 
Therefore, once such timescale is calibrated through 
the observed ratio SNRi a (2 = 0)/SNRi a (2 = 1), one 


should be able to constrain the star formation his¬ 
tory of the early universe by comparing the predicted 
SN la rate at 2 > 2 with the observations. 

The main results of this exercise can be summarized 
as follows: 

• At the present epoch, the predicted Type 
II(+Ib/c) frequency matches remarkably well 
the observed local value. Note, however, that 
rates obtained from traditional distant (beyond 
4 Mpc) sample might need to be increased by 
a factor of 1.5-2 because of severe selection ef¬ 
fects against Type II’s fainter than My = —16 
(Woltjer 1997).' 

• In the interval 0 < 2 < 1, the predicted rate of 
SN la is a sensitive function of the characteristic 
delay timescale between the collapse of the pri¬ 
mary star to a WD and the SN event. Accurate 
measurements of SN rates in this redshift range 
will improve our understanding of the nature of 
SN la progenitors and the physics of the explo¬ 
sions. Ongoing searches and studies of distant 
SNe should soon provide these rates, allowing a 
universal calibration of the Type la phenomenon. 

• While Type la rates at 1 < 2 < 2 will offer valu¬ 
able information on the star formation history 
of the universe at earlier epoch, the full picture 
will only be obtained with statistics on Type la 
and II SNe at redshifts 2 < z < 4 or higher. 
At these epochs, the detection of Type II events 
must await NGST. A SN II has a typical peak 
magnitude Mb « —17: placed at 2 = 3, such an 
explosion would give rise to an observed flux of 
15 nJy (assuming a flat cosmology with qo = 0.5 
and Hq = 50 h .50 kms -1 Mpc -1 ) at 1.8 /im. At 
this wavelength, the imaging sensitivity of an 8m 
NGST is 1 nJy (10 4 s exposure and 10er detec¬ 
tion threshold), while the moderate resolution 
(A/AA = 1000) spectroscopic limit is about 50 
times higher (10 s s exposure per resolution el¬ 
ement and 10 a detection threshold) (Stockman 
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et al. 1998). The several weeks period of peak 
rest-frame blue luminosity would be stretched 
by a factor of (1 + z) to few months. Figure 3 
shows the cumulative number of Type II events 
expected per year per 4' x 4' field. Depending on 
the history of star formation at high redshifts, 
NGST should detect between 7 (in the merging 
model) and 15 (in the monolithic collapse sce¬ 
nario) Type II SNe per field per year in the in¬ 
terval 2 < 2 < 4. The possibility of detecting 
Type II SNe at z > 5 from an early popula¬ 
tion of galaxies has been investigated by Miralda- 
Escude & Rees (1997). By assuming these are 
responsible for the generation of all the metals 
observed in the Lya forest at high redshifts, a 
high baryon density (Df,ft| 0 = 0.1), and an av¬ 
erage metallicity of 0.01 Z Q , Miralda-Escude & 
Rees estimate NGST should observe about 16 SN 
II per field per year with z > 5. Note, however, 
that a metallicity smaller by a factor ~ 10 com¬ 
pared to the value adopted by these authors has 
been recently derived by Songaila (1997). For 
comparison, the models discussed above predict 
between 1 and 10 Type II SNe per field per year 
with z > 4. 


3. WHAT KEEPS THE UNIVERSE IONIZED 
AT Z = 5? 


The existence of a filamentary, low-density inter- 
galactic medium (IGM), which contains the bulk of 
the hydrogen and helium in the universe, is predicted 
as a product of primordial nucleosynthesis and of 
hierarchical models of gravitational instability with 
“cold dark matter” (CDM) (Gen et al. 1994; Zhang 
et al. 1995; Hernquist et al. 1996). The application 
of the Gunn-Peterson constraint on the amount of 
smoothly distributed neutral material along the line 
of sight to distant objects requires the hydrogen com¬ 
ponent of the diffuse IGM to have been highly ion¬ 
ized by z ~ 5 (Schneider et al. 1991), and the he¬ 
lium component by z « 2.5 (Davidsen et al. 1996). 
The plethora of discrete absorption systems which 
give origin to the Lya forest in the spectra of back¬ 
ground quasars are also inferred to be strongly pho- 
toionized: we know from QSO absorption studies, in 
fact, that neutral hydrogen accounts for only a small 
fraction, ~ 10%, of the nucleosynthetic baryons at 
early epochs (Lanzetta et al. 1995). It thus appears 
that substantial sources of ultraviolet photons were 
present at z > 5, perhaps low-luminosity quasars 
(Haiman & Loeb 1998) or a first generation of stars in 
virialized dark matter halos with T v j r ~ 10 4 — 10 5 K 
(Couchman & Rees 1986; Ostriker & Gnedin 1996; 
Haiman & Loeb 1997; Miralda-Escude & Rees 1998): 
early star formation provides a possible explanation 
for the widespread existence of heavy elements in the 
IGM (Cowie et al. 1995). More in general, establish¬ 
ing the character of cosmological ionizing sources is 
an efficient way to constrain competing models for 
structure formation in the universe, and to study the 
collapse and cooling of small mass objects at early 
epochs. 



Figure 4- Comoving space density of bright QSOs as a 
function of redshift. The data points with error bars are 
taken from HS (filled dots), WHO (filled squares), SSG 
(crosses), and KDC (filled pentagon). The points have 
been normalized to the z = 2.5 space density of quasars 
with Mb < —26 (Mb < —27 in the case of KDC). The 
empty triangles show the space density (normalized to the 
peak) of the Parkes flat-spectrum radio-loud quasars with 
P > 7.2 x 10 26 W Hz- 1 sr” 1 (Hook et al. 1998). 


3.1. Quasars 


A decline in the space density of bright quasars at 
redshifts beyond ~ 3 was first suggested by Osmer 
(1982), and has been since then the subject of a 
long-standing debate. In recent years, several optical 
surveys have consistently provided new evidence for 
a turnover in the QSO counts (Hartwick & Schade 
1990, HS; Warren et al. 1994, WHO; Schmidt et al. 
1995, SSG; Kennefick et al. 1995, KDC). The in¬ 
terpretation of the drop-off observed in optically se¬ 
lected samples is equivocal, however, because of the 
possible bias introduced by dust obscuration arising 
from intervening systems (Ostriker & Heisler 1984). 
Radio emission, on the other hand, is unaffected by 
dust, and it has recently been shown (Shaver et al. 
1996) that the space density of radio-loud quasars 
also decreases strongly for z > 3 (Fig. 4), demon¬ 
strating that the turnover is indeed real and that 
dust along the line of sight has a minimal effect on 
optically-selected QSOs. The QSO emission rate of 
hydrogen ionizing photons per unit comoving volume, 
A f, has been recently recalculated by Madau et al. 
(1998b), and is shown in Figure 5. It is important 
to notice that the procedure adopted to derive this 
quantity implies a large correction for incomplete¬ 
ness at high-z. With a fit to the quasar luminosity 
function (LF) which goes as <t>(L) oc L~ 164 at the 
faint end (Pei 1995), the contribution to the emissiv- 
ity converges rather slowly, as L 0 36 . At z = 4, for 
example, the blue magnitude at the break of the LF 
is M* ss —25.4, comparable or slightly fainter than 
the limits of current high-z QSO surveys. A large 
fraction, about 90% at z = 4 and even higher at ear¬ 
lier epochs, of the ionizing emissivity in our model is 
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therefore produced by quasars that have not been ac¬ 
tually observed, and are assumed to be present based 
on an extrapolation from lower redshifts. The value 
of Af obtained by including the contribution from ob¬ 
served quasars only would be much smaller at high 
redshifts than shown in Figure 5. 


3.2. Star-forming Galaxies 

Galaxies with ongoing star-formation are another ob¬ 
vious source of Lyman continuum photons. A com¬ 
posite UV luminosity function of Lyman-break galax¬ 
ies at z ~ 3 has been recently derived by Dick¬ 
inson (1998). It is based on spectroscopically and 
photometrically selected galaxies from the ground- 
based and HDF samples, and spans about a fac¬ 
tor 50 in luminosity from the faint to the bright 
end. Because of the uncertanties that still remain 
in the rescaling of the HDF data points to the 
ground-based data, the comoving luminosity density 
at 1500 A is estimated to vary within the range 1.6 to 
3.5 x 10 26 ergss -1 Hz -1 Mpc -3 . Since the rest-frame 
UV continuum at 1500 A (redshifted into the visi¬ 
ble band for a source at z ~ 3) is dominated by the 
same short-lived, massive stars which are responsible 
for the emission of photons shortward of the Lyman 
edge, the needed conversion factor, about one ion¬ 
izing photon every 10 photons at 1500 A, is fairly 
insensitive to the assumed IMF and is independent 
of the galaxy history for t 10 73 yr. 

Figure 6 shows the Lyman continuum luminosity 
function of galaxies at z « 3 [at all ages >0.1 Gyr one 
has L(1500)/L(912) ss 6 for a Salpeter mass function 
and constant star formation rate], compared to the 
distribution of QSO luminosities at the same redshift. 
The comoving ionizing emissivity due to Lyman- 
break galaxies is 4.2± 1.5 x 10 25 ergs s -1 Hz -1 Mpc -3 , 
between 2 and 4 times higher than the estimated 
quasar contribution at z = 3. This number neglects 
any correction for intrinsic H i absorption: the data 
points plotted in Figure 6 assume a value of / esc = 0.5 
for the unknown fraction of ionizing photons which 
escapes the galaxy H i layers into the intergalactic 
medium. 

The LF of Lyman-break galaxies at 0 > 4 is highly 
uncertain. An analysis of the H-band dropouts in 
the HDF - candidate star-forming objects at 3.5 < 
z < 4.5 - seems to imply a decrease in the comov¬ 
ing UV galaxy emissivity by about a factor of 2.5 
in the interval 2.75 < z < 4 (Madau et al. 1996, 
1998c). In this sense star-forming galaxies with SFR 
in excess of 0.5MQyr -1 may have a negative evolu¬ 
tion with lookback time similar to the one observed 
in bright QSOs, but the error bars are still rather 
large. Adopting a L(1500) to L(912) conversion fac¬ 
tor of 6, we estimate a comoving ionizing emissivity 
of 1.7 ± 1.1 X 10 25 / esc ergs s -1 Hz -1 Mpc -3 at z « 4. 
One should note that, while highly reddened sub-mm 
sources at high redshifts (Lilly, this volume) would 
be missed by the dropout color technique (which iso¬ 
lates sources that have blue colors in the optical and 
a sharp drop in the rest-frame UV), it seems un¬ 
likely that very dusty objects (with / esc <C 1) would 
contribute in any significant manner to the ionizing 
metagalactic flux. 



Figure 5. Comoving emission rate of hydrogen Lyman 
continuum photons (solid line) from QSOs, compared with 
the minimum rate (dashed line) which is needed to fully 
ionize a fast recombining (with clumping factor C = 30) 
universe with D{,/i 2 0 = 0.08. Models based on photoion¬ 
ization by quasar sources appear to fall short at z = 5. 
The data points with error bars show the estimated con¬ 
tribution of star-forming galaxies at Z ~ 3 and , with sig¬ 
nificantly larger uncertainties, at z « 4. The fraction 
of Lyman continuum photons which escapes the galaxy 
H I layers into the intergalactic medium is taken to be 
/esc = 0.5. 


3.3. Reionization of the Universe 


In inhomogeneous reionization scenarios, the his¬ 
tory of the transition from a neutral IGM to one 
that is almost fully ionized can be statistically de¬ 
scribed by the evolution with redshift of the vol¬ 
ume filling factor or porosity Q(z) of H II , He II , 
and He in regions. The radiation emitted by spa¬ 
tially clustered stellar-like and quasar-like sources - 
the number densities and luminosities of which may 
change rapidly as a function of redshift - coupled 
with absorption processes in a medium with a time- 
varying clumping factor, all determine the complex 
topology of neutral and ionized zones in the universe. 
When Q « 1 and the radiation sources are ran¬ 
domly distributed, the ionized regions are spatially 
isolated, every UV photon is absorbed somewhere 
in the IGM, and the ionization process cannot be 
described as due to a statistically homogeneous ra¬ 
diation field. As Q grows, the crossing of ionization 
fronts becomes more and more common, and the neu¬ 
tral phase shrinks in size until the reionization pro¬ 
cess is completed at the “overlap” epoch, when every 
point in space is exposed to Lyman continuum radi¬ 
ation. The filling factor of H n regions in the uni¬ 
verse, Qhii, is equal at any given instant t to the 
integral over cosmic time of the rate of ionizing pho¬ 
tons emitted per hydrogen atom and unit cosmologi¬ 
cal volume by all radiation sources present at earlier 
epochs, f* fii on (t')dt'/nn(t'), minus the rate of radia¬ 
tive recombinations, f* Qmi(t')dt'/i Tec (t') (Madau et 











7 



M #b (912A) 


Figure 6. The 912 ^4 luminosity function of galaxies at 
z ~ 3 (solid line), compared to the distribution of QSO 
luminosities at the same redshift (dashed line). The for¬ 
mer assumes a Salpeter IMF with constant constant star 
formation rate (age=l Gyr): -Mab(912A)= —19 corre¬ 
sponds to a rate of 13M0yr _1 . The solid and dashed 
lines represent functional fits to the data points , and the 
dotted lines their extrapolation. 


a definite picture for the origin of intervening absorp¬ 
tion systems, one of an interconnected network of 
sheets and filaments, with virialized systems located 
at their points of intersection. In the simulations of 
Gnedin & Ostriker (1997), for example, the clumping 
factor rises above unity when the collapsed fraction 
of baryons becomes non negligible, i.e. z < 20, and 
grows to C > 10 (40) at 2 ~ 8 (5) (because of fi¬ 
nite resolution effects, numerical simulations will ac¬ 
tually underestimate clumping): the recombination 
timescale is much shorter than that for a uniform 
IGM, and always shorter than the expansion time. 

For an IGM with C > 1, and in the case nj on and 
C do not vary rapidly over a timescale t rec , one can 
expand around t to find 

Qmi(t) « ^f- rec . (7) 

n H 

The porosity of ionized bubbles is then approximately 
given by the number of ionizing photons emitted per 
hydrogen atom in one recombination time. In other 
words, because of hydrogen recombinations, only a 
fraction t lec /t (~ a few per cent at z = 5) of the 
photons emitted above 1 ryd is actually used to ion¬ 
ize new IGM material. The universe is completely 
reionized when Qhii = 1, i.e. when 

h'ionGec = tl [ ]. (8) 


al. 1998b). Differentiating one gets 

dQmi _ fhon Qhii ^ 

dt 7T-H tree 

It is this simple differential equation - and its equiv¬ 
alent for expanding helium zones - that statistically 
describes the transition from a neutral universe to 
a fully ionized one, independently, for a given emis- 
sivity, of the complex and possibly short-lived emis¬ 
sion histories of individual radiation sources, e.g., on 
whether their comoving space density is constant or 
actually varies with cosmic time. Here hi on is the 
emission rate of ionizing photons per unit proper vol¬ 
ume, n h is the mean hydrogen density of the expand¬ 
ing IGM, nn(O) = 1.7 x 10“' (Hf,/i§ 0 /0.08) cm -3 , 

t rec = [(1+2 X )n H Q B C 1 ] -1 = 0.3 Gyr C^ 1 

( 6 ) 

is the volume-averaged gas recombination timescale, 
as is the recombination coefficient to the excited 
states of hydrogen, \ the helium to hydrogen cos¬ 
mic abundance ratio, C = {'«-hii)/^-hii is the ionized 
hydrogen clumping factor, and a gas temperature of 
10 4 K has been assumed. Clumps which are dense 
and thick enough to be self-shielded from UV radia¬ 
tion will stay neutral and will not contribute to the 
recombination rate. 

An empirical determination of the dumpiness of the 
IGM at high redshifts is hampered by our poor 
knowledge of the ionizing background intensity and 
the typical size and geometry of the absorbers. Nu¬ 
merical N-body/hydrodynamics simulations of struc¬ 
ture formation in the IGM within the framework of 
CDM dominated cosmologies have recently provided 


4. A POPULATION OF EARLY DWARFS 


As tree Cf at high redshifts, it is possible to com¬ 
pute at any given epoch a critical value for the photon 
emission rate per unit cosmological comoving volume, 
J\f c . independently of the (unknown) previous emis¬ 
sion history of the universe: only rates above this 
value will provide enough UV photons to ionize the 
IGM by that epoch. One can then compare our deter¬ 
minations of N c to the estimated contribution from 
QSOs and star-forming galaxies. Equation (|j) can 
then be rewritten as 


tre c{z) 


(IQ 512 s _1 Mpc _3 )C 3 o 



(9) 

The uncertainty on this critical rate is difficult to es¬ 
timate, as it depends on the clumping factor of the 
IGM (scaled in the expression above to the value 
inferred at z = 5 from numerical simulations) and 
the nucleosynthesis constrained baryon density. A 
quick exploration of the available parameter space 


indicates that the uncertainty on J\f c could easily be 
of order ±0.2 in the log. The evolution of the critical 
rate as a function of redshift is plotted in Figure 5. 
While I\f c is comparable to the quasar contribution at 
2 > 3, there is some indication of a significant deficit 
of Lyman continuum photons at z = 5. For bright, 
massive galaxies to produce enough UV radiation at 
2 = 5, their space density would have to be compara¬ 
ble to the one observed at 2 « 3, with most ionizing 
photons being able to escape freely from the regions 
of star formation into the IGM. This scenario may 
be in conflict with direct observations of local star- 
bursts below the Lyman limit showing that at most a 








few percent of the stellar ionizing radiation produced 
by these luminous sources actually escapes into the 
IGM (Leitherer et al. 1995). If, on the other hand, 
faint QSOs with (say) Mab = —19 at rest-frame ul¬ 
traviolet frequencies were to provide all the required 
ionizing flux, their comoving space density would be 
such (0.0015 Mpc -3 ) that about 50 of them would 
expected in the HDF down to Iab = 27.2. At z > 5, 
they would appear very red in V — I as the Lya forest 
is shifted into the visible. This simple model can be 
ruled out, however, as there is only a handful (7) of 
sources in the HDF with (V — I) ab > 1-5 mag down 
to this magnitude limit. 


It is interesting to convert the derived value of J\f c 
into a “minimum” star formation rate per unit (co¬ 
moving) volume, p*): 

P*( z ) = tic{z) x 10 -53 ' 1 / -1 (10) 

« 0.013/ e - ! (-y-) M 0 yr -1 Mpc -3 (11) 

(Salpeter IMF with solar metallicity), comparable 
with the value directly “observed” (i.e., uncorrected 
for dust reddening) at z ~ 3 (Dickinson 1998). The 
same massive stars that dominate the Lyman contin¬ 
uum flux also manufacture and return most of the 
metals to the ISM. In the approximation of instanta¬ 
neous recycling, the rate of ejection of newly sinthe- 
sized heavy elements which is required to keep the 
universe ionized at redshift z is, from equation ©> 

pz{z) = y(l- R)p*{z) >3.5 x 1(T 4 (12) 

^ lM0 V- lM PC- a . (U) 

where y is the net, IMF-averaged “yield” of returned 
metals, Zq = 0.02, and R ~ 0.3 is the mass frac¬ 
tion of a generation of stars that is returned to the 
interstellar medium. At z = 5, and over a timescale 
of At = 0.5 Gyr (corresponding to a formation red- 
shift Zf = 10) such a rate would generate a mean 
metallicity per baryon in the universe of 


(Z) 


8nGpz(5)At 
3 H*Sl b 


> 0.002 



/esc Zq ; (14) 


comparable with the level of enrichment observed 
in the Lya forest at z ~ 3 (Songaila 1997): more 
than 2% of the present-day stars would need to have 
formed by z ~ 5. It has been recently suggested 
(Miralda-Escude & Rees 1998) that a large num¬ 
ber of low-mass galactic halos, expected to form at 
early times in hierarchical clustering models, might 
be responsible for photoionizing the IGM at these 
epochs. According to spherically-symmetric simu¬ 
lations (Thoul et al. 1996), photoionization heating 
by the UV background flux that builds up after the 
overlapping epoch completely suppresses the cooling 
and collapse of gas inside the shallow potential wells 
of halos with circular velocities < 35 kms -1 . Halos 
with circular speed u c i rc = 50 kms -1 , correspond¬ 
ing in top-hat spherical collapse to a virial temper¬ 
ature T v j r = 0.5 pm p v1njk « 10 5 K and halo mass 

M = 0.1i4 c /GH « 4 x 10 9 [(1 + z)/6\~ 3 / 2 h^ M 0 , 
appear instead largely immune to this external feed¬ 
back (but see Navarro & Steinmetz 1997). In these 


systems rapid cooling by atomic hydrogen can then 
take place and a significant fraction, /fit,, of their to¬ 
tal mass may be converted into stars over a timescale 
At comparable to the Hubble time (with / close 
to 1 if the efficiency of forming stars is high). If 
high- 2 : dwarfs with star formation rates ffl^M/At ~ 
0.3/o. 5 A^ 5 M 0 yr -1 were actually responsible for 
keeping the universe ionized at 2 ~ 5, their comoving 
space density would have to be 

0.013/ -1 M 0 yr -1 Mpc -3 / /esc/ y 1 M 3 

0.6/ M 0 yr -1 V 0-25 J P ’ 

(15) 

two hundred times larger than the space den¬ 
sity of present-day galaxies brighter than L*(4400), 
and about five hundred times larger than that 
of Lyman-break objects at 2 k 3 with M < 
M£ b ( 1500), i.e. with star formation rates in excess 
of 10 M 0 yr -1 . Only a rather steep luminosity func¬ 
tion, with Schechter slope a ~ 2, would be consistent 
with such a large space density of faint dwarfs and, 
at the same time, with the paucity of brighter B- and 
V-band dropouts observed in the HDF. The number 
density on the sky would be « 0.2 arcsec -2 , corre¬ 
sponding to more than three thousands sources in the 
HDF. With a typical apparent magnitude at z = 5 
of /ab ~ 29.5 mag (assuming / = 0.5), these are too 
faint to be detected by the HST, but within the range 
of the proposed NGST. A higher density of sources 
- which would therefore have to originate from lower 
amplitude peaks - would be required if the typical ef¬ 
ficiency of star formation and/or the escape fraction 
of ionizing photons were low, (/,/ esc ) « 1- In this 
case the dwarfs could still be detectable if a small 
fraction of the gas turned into stars in very short 
bursts. 
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